Electronic transport and quantum Hall effect in bipolar graphene p-n-p junctions 
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We have developed a device fabrication process to pattern graphene into nanostructures of arbi- 
trary shape and control their electronic properties using local electrostatic gates. Electronic trans- 
port measurements have been used to characterize locally gated bipolar graphene p-n-p junctions. 
We observe a series of fractional quantum Hall conductance plateaus at high magnetic fields as the 
local charge density is varied in the p and n regions. These fractional plateaus, originating from 
chiral edge states equilibration at the p-n interfaces, exhibit sensitivity to inter-edge backscattering 
which is found to be strong for some of the plateuas and much weaker for other plateaus. We use this 
effect to explore the role of backscattering and estimate disorder strength in our graphene devices. 



O 

(N 

> 

O 

cn 

in 
o 



Graphene, a recently discovered single sheet of 
graphite [1], stands out as an exceptional material both 
in terms of the fundamental physics associated with its 
unique " quasi-relativistic" carrier dynamics and poten- 
tial applications in electronic devices Most exper- 
iments to date depended on the presence of a heav- 
ily doped Si substrate which serves as a global back 
gate, inducing charge density via the electric field ef- 
fect. Although such global gate approach yields inter- 
esting transport phenomena [H, 0, 13, H, 0, Q , it repre- 
sents only the first step towards more complex graphene 
devices. The use of local gates enables the fabrication 
of in-plane graphene heterostructures [1, [3]. A num- 
ber of applications for local gate devices have been pro- 
posed in order to investigate Klein tunneli ng p , electron 
Veselago-lens |10], quantum point contacts [ll^ and quan- 
tum dots [12, 13]. 

In this letter, we study locally gated graphene devices 
in the quantum Hall (QH) regime. By independently 
varying voltage on the back gate and local gate, we can 
study bipolar QH transport in graphene p-n-p hetero- 
j unctions in different charge density regimes. We find 
a series of fractional QH conductance plateaus as the 
local charge density is varied in the p and n regions. 
Recently similar QH effect has been reported in a sin- 
gle p-n graphene hetero junction [15]. Our double junc- 
tion system allows to study new interesting transport 
regimes that are absent in the QH edge transport in 
a single junction, in particular, partial equilibration of 
graphene QH edge states. Conspicuously, some of our 
fractionally quantized plateaus are found to be consider- 
ably more fragile with respect to disorder than the other 
plateaus. We analyze the distinction in roughness of dif- 
ferent plateaus and show that it points to the importance 
of inter-edge backscattering in our narrow graphene sam- 
ples. 

Our devices have been fabricated following a com- 
bination of an etching and dielectric deposition tech- 
niques [l^. Fig. 1(a) shows examples of nanostructured 
graphene devices after the complete fabrication process. 




FIG. 1: (a) Scanning electron microscopy (SEM) picture 
showing several complete two-probe devices with local gates. 
Scale bar 2 /xm. The electrodes contact the widest parts of 
the structure. The local gates cover the graphene constric- 
tions and part of the central graphene bars, (b) Schematic 
side view of the devices. Graphene is contacted by source (S) 
and drain (D) electrodes, and separated from the back gate 
plane by 300 nm Si02, and from the local gate (LG) by the 
top dielectric. The back gate (with voltage Vbg) is coupled 
to the entire graphene structure via the capacitor Cbg- The 
local gate only couples to part of the structure via Clg- 



The devices consist of a designed graphene nanostructure 
sandwiched between two dielectrics, with a global back 
gate (the highly doped Si substrate) and partially covered 
by one or more local gates (Fig. 1(b)). Such gates allow 
us to tune the location of the Fermi energy in graphene 
globally (via the back gate voltage, Vbg) or locally (via 
the local gate voltage, Vlg)- The conductance, G, of 
our devices is measured at cryogenic temperatures (1.5 
to 4.2 K), as a function of Vbg and Vlg, by using a lock- 
in technique with an ac excitation voltage of 100 jiY. 
The single layer character of our devices is determined 
by Raman spectroscopy [17] and/or quantum Hall effect 
measurements [3, 4]. 

Bulk graphene is a zero band gap semiconductor. 
Therefore, the Fermi energy in graphene can be contin- 
uously varied from valence to conduction band via the 
field effect. Incorporating local gates allows us to induce 
different charge densities at different sample regions. Of 
particular interest is the case when the Fermi energy in 
one region is in the valence band (p-type) while in the 
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FIG. 2: (a) G{Vbg) for a graphene p-n-p junction, extractec 
from (b), showing the two conductance minima associatec 
with Dirac valleys in the graphene leads and under the lo- 
cal gate. Inset: false color SEM picture of a patternec 
graphene bar with contacts and local gate. Scale bar rep- 
resents 1 /xm. (b) Two-dimensional plot of G{Vlg,Vbg) foi 
the device shown in (a). 



other region it is in the conduction band (n-type). 

We have fabricated graphene p-n-p/n-p-n devices with 
different gated channel width and length (see Fig. 1(a)). 
In total 6 devices in two different single layer graphene 
pieces were studied and found to exhibit similar char- 
acteristics. Fig. 2(b) shows G(Vlg, Vbg) at zero mag- 
netic field for a typical device [14]. The most prominent 
feature is the presence of two conductance minima val- 
leys: one approximately horizontal, independent of the 
voltage Vlg^ and the other diagonal. The first valley 
tracks the charge neutrality point (or, Dirac point) in 
the regions outside the local gate, further denoted as 
graphene leads (GLs).The second valley, whose position 
is controlled both by the local and the back gate volt- 
age, tracks the neutrality point in the region under the 
local gate (LGR). The slope of this valley is equal to the 
ratio of the local capacitance to the back gate capaci- 
tance. For the device shown in Fig. 2(b) this ratio is 
10.5. Note that the typical local gate breakdown volt- 
age for our dielectric heterostructure is larger than 12 V, 
which, together with a ten-fold enhancement of the gate 
coupling factor, means that we can induce local charge 
densities at least comparable, and often larger, than with 
the back gate (where leakage typically starts to occur at 
Vbg - 100 V). 

Away from the Dirac valleys, the conductance in- 
creases with increasing charge density (see Fig. 2). The 
two valleys separate four regions in the plot: p-p'-p^ 
p-n'-p^ n-n'-n and n-p'-n^ where n/p refers to nega- 
tive/positive charge density and the prime indicates den- 
sity in the LGR. The conductance is not symmetric across 
the valleys, because for opposite polarities there is an ex- 
tra contribution due to the resistance of the two p-n in- 
terfaces [14] . We also note here that even in the p-n-p and 
n-p-n regions, the device shows considerable conductance 
(G > /h) without any signature of rectifying behavior, 
as expected for transport in a zero-gap hetero junction. In 
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FIG. 3: (a) Color map of conductance GiyLC^VBo) at mag- 
netic field 5 = 13 T, and T = 4.2 K. The black cross indi- 
cates the location of filling factor zero in LGR and GLs. Inset: 
Conductance at zero B in the same (Vlg, Vbg) range and the 
same color scale as main figure (white denotes G > 10.5e^//i). 
(b) G{Vlg) extracted from (a), red trace, showing fractional 
values of the conductance. Numbers on the right indicate 
expected fractions for the various filling factors (red num- 
bers indicate the filling factor, iy\ in LGR), see also (g). (c) 
G{Vlg) (projection of orange trace from (a) onto VLG-axis). 
Orange numbers indicate filling factor, in the GLs. (d) to 
(f): different edge state diagrams representing possible equi- 
libration processes taking place at different charge densities 
in the GLs and LGR. The purple region indicates the LGR. 
Yellow boxes indicate contact electrodes, (g) Simulated color 
map of the theoretical conductance plateaus expected from 
the mechanisms shown in (d-f) for different filling factors in 
the GLs and LGR. The numbers in the rhombi indicate the 
conductance at that plateau. The color scale is identical to 
that of (a). 



fact, graphene is the only two-dimensional electron gas 
(2DEG) where in-plane bipolar heterostructures p-n'-p 
and n-p'-n can be studied in the linear response regime. 

The lateral graphene hetero junctions exhibit interest- 
ing phenomena at high magnetic field. One of the hall- 
marks of graphene is the relativistic integer Quantum 
Hall (QH) effect, manifested in a series of conductance 
plateaus at half-integer multiples of 4e'^/h [1,0, Such 
unique QH plateau structure can be attributed to an odd 
number of QH edge states that carry current with con- 
ductance 2e^/h [2C 
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The capability of placing lo- 
cal tunable electrostatic barriers/wells along the current 
pathway allows us to use QH mode propagation to ex- 
plore intrinsic transport characteristics of graphene het- 
eroj unction structures. 

Fig|3^ shows G{Vlg^Vbg) at high magnetic field, 
for a typical device. In addition to the four regions 
seen at zero field (FiglSfci, inset), the conductance map 
reveals a rich pattern consisting of rhombi and bands 
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where G exhibits plateaus. Over ah this pattern is sym- 
metric with respect to the neutrahty point (marked with 
a black cross) which corresponds io v = v' = {). (Here 
V and Tj' are the Landau level filling factors, equal to 
Tice/hB with Uc the carrier density in GLs and LGR, re- 
spectively.) While the plateaus in conductance at high 
densities > 6, are well accounted for by a two 

resistors-in-series model, with each resistor correspond- 
ing to the QH conductance of GLs, a more complex and 
interesting behavior is observed at lower densities, where 
resistors cease to add up in a classical fashion. 

The nonclassical behavior is found in particular at low 
filling factors, especially when either or u equals +2 
or —2 (red and orange traces in Fig|3l(a)). Notably, we 
observe conductance plateaus at values close to fractional 
values of the conductance quantum, e'^/h. Such simple 
fractions include, for example, (2/3)e^//i, (6/7)e^//i and 
(10/9)e^//i (FiglSJb)). These values are in sharp contrast 
to the conductance plateaus at (2,6, 10...)e^//i, observed 
in homogenous two-terminal devices [7|. 

The unusual fractional conductance plateau patterns 
can be analyzed by using models developed for QHE 
mode propagation in 2DEGs with density gradients [22, 
[isj . Our graphene system, however, represents a distinct 
advantage owing to the possibility of creating opposite 
polarities of charge carriers in adjacent regions. 

The simplest case to consider is when the polarity of 
GLs and LGR is the same (either n-n'-n or p-p'-p), and 
the LGR density is lower than the GLs density: < 
In this case, as shown in Fig|3l(d), the number of QH 
edge modes is larger in the GLs than in the LGR. The 
modes existing only in the GLs are fully reflected at the 
LGR-GLs interfaces [23| , while those present in both 
regions exhibit full transmission, giving rise to the net 
conductance G = {e'^/h)\u^\. 

A more interesting situation occurs when the LGR den- 
sity is higher than the GLs density, with the LGR and 
GLs polarities still the same. In this case the number of 
edge states is smaller in the GLs than in the LGR (see 
FigEl^e)). Crucially, the states circulating in LGR can 
produce partial equilibration among the different chan- 
nels, because they couple modes with different electro- 
chemical potentials. To analyze this regime, we suppose 
that current / is injected from the left lead, while no 
current is injected from the right lead. Then the conser- 
vation of current yields I -\- 14 = /i, h = h (the LGR 
edges are labeled by 1, 2, 3, 4 as shown on Figl3l(e)). 
Assuming that the current at the upper and lower LGR 
edges is partitioned equally among available edge modes, 
we obtain the relations for the current flowing out of 
these edges: I2 = r/i, I4 = r/3, (r = 1 — I'/i'')- Solving 
these equations for /1...4, we determine the current flow- 
ing in the drain lead as lout = h — h and flnd the net 
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(l^'l > IH), (1) 



where v^v' = ±2, ±6.... We emphasize that this partial 
equilibration regime can only occur in the presence of two 
n-n' or p-p' interfaces, and would not occur in a single 
n-n' or p-p' junction [15I. [18] . 

The last, but most unique case is when the GLs and 
LGR have opposite carrier polarity. In this case, the edge 
states counter-circulate in the p and n areas, running par- 
allel to each other along the p-n interface (see FiglSl^f)). 
Such propagation, leading to mixing among edge states, 
results in full equilibration at the p-n interfaces:/i = r/2, 
Is = r/4, (r = |i^'|/(|i^| + Combining this with 

current conservation, in this case written as / + /i = /4, 
I2 = ^3, we flnd the currents and obtain the conductance 
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{uu' < 0), (2) 



where v' = ±2, ±6.... The net conductance in this case 
is described by three quantum resistors in series. 

The summary of all possible conductance values for 
these three regimes is shown as a color map in Figl3^. 
Our flrst observation is that the structure of the ex- 
perimental pattern resembles qualitatively the theoret- 
ical one when the fllling factor equals ±2 either in the 
GLs or in the LGR. For a quantitative analysis, we 
choose two cuts extracted from Figl3fa), showing con- 
ductance for flxed u = -2 (FigEb) and = 2 (FigEt). 
We record reasonably good plateaus at G = (2/3)e^//i, 
G = (10/9)e^//i as well as other fractions discussed 
above, with the only exception of a considerably more 
poor plateau with G = 2e^/h (see below). Of partic- 
ular interest is the non-monotonic conductance behav- 
ior in Figlajb) for v' = 2,-2,-6,-10 (with u = -2), 
which reflects the full equilibration edge state trans- 
mission partial equilibration sequence. This is in con- 
trast with the monotonic behavior of G in Fig|3l^c) for 
u = —2 ^ 2 ^ 6 (with = 2), where only the full 
equilibration and full transmission regimes are expected. 

We have measured in total four devices which all ex- 
hibit similar conductance patterns. Notably, in none of 
these devices G reached the full 2e^//i value at z^', = ±2, 
whereas other conductance plateaus were well devel- 
oped. The lack of quantization points to the presence 
of backscattering between opposite edges of our sample, 
which may occur in LGR bulk or in the transitional re- 
gions at the LGR-GLs junctions. 

Why are the u'^u = ±2 plateaus so sensitive to 
backscattering, while other plateaus are not? To gain 
insight into this question, we now investigate how robust 
are the results ([I]),© with respect to bulk conduction in 
our QH system. To that end, we consider a 2D transport 
model describing the system by local conductivity. Here 
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FIG. 4: Two-terminal conductance of a lateral heteroj unction 
in a QH state (see inset) as a function of the longitudinal 
conductivity under local gate (LGR) for the states with v — 
—2 and v' — ±2,±6, — 10 (red trace in Fiffj3la,b)). Results 
obtained from the 2D transport model [2^, 12511261] for LGR of 
size 500 nm x 700 nm are shown. Finite Uxx has considerable 
effect on the conductance of the state with v' — —2, but 
very little effect on other states. This explains the difference 
in roughness of the observed plateaus for v' ^ —2 and v' — 
-2. 

we focus on the simplest situation, taking the longitudi- 
nal conductivity Gxx nonzero in the gated region (LGR) 
and zero outside (GLs), and the Hall conductivity Oxy 
equal to v'e^jh {ye^ jh) for LGR (GLs). 

An exact solution for 2D current and potential distri- 
bution for this problem was obtained [24] by adapting the 
conformal mapping technique developed in Refs.[25, 26]. 
The resulting two-terminal conductance G of the frac- 
tionally quantized states from the v = —2 trace (Figl3)3) 
are displayed in FiglH First, we note that the limit- 
ing values of G at g^x agree with the simple frac- 
tions ([I]),© derived above. Furthermore, the effect of 
finite a^x is considerably stronger for the v = v' = —2 
state than for all other states — it is linear rather than 
quadratic at small g^x- Comparing to the deviation from 
the quantized value in FiglSl we estimate g^x ^ 0.5e^//i. 

We thus infer that weak backscattering is non- 
detrimental for all the states except v = v' ^ which is 
in agreement with the observed stability of fractional 
plateaus. This conclusion also agrees with the general in- 
tuition that current paths in a QH system are constrained 
stronger when density is varying in space than when it is 
constant [27]. We therefore believe that our understand- 
ing of stability of the observed fractional plateaus is quite 
generic and insensitive to whether the backscattering in 
our graphene devices occurs mainly in LGR bulk or at 
the LGR-GLs interfaces. 
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